Test for anisotropy in the mean of the CMB temperature fluctuation in spherical 
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The standard models of inflation predict statistically homogeneous and isotropic primordial fluc- 
tuations, which should be tested by observations. In this paper we illustrate a method to test the 
statistical isotropy of the mean of the cosmic microwave background temperature fluctuations in 
the spherical harmonic space and apply the method to the Wilkinson Microwave Anisotropy Probe 
seven-year observation data. A classical method to test a mean, like the simple Student's t test, is 
not appropriate for this purpose because the Wilkinson Microwave Anisotropy Probe data contain 
anisotropic instrumental noise and suffer from the effect of the mask for the foreground emissions 
which breaks the statistical independence. Here we perform a band-power analysis with Monte Carlo 
simulations in which we take into account the anisotropic noise and the mask. We find evidence 
of a non-zero mean at 99.93 % confidence level in a particular range of multipoles. The evidence 
against the zero-mean assumption as a whole is still significant at the 99 % confidence level even if 
the fact is taken into account that we have tested multiple ranges. 

PACS numbers: 98.70.Vc, 98.80.-k, 98.80.Es 
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I. INTRODUCTION 

Inflation provides a successful mechanism of generat- 
ing primordial density perturbations that give rise to 
the large-scale structure of the Universe and temperature 
anisotropy in the cosmic microwave background (CMB). 
One of the important consequences of inflation models, 
which should be tested with observations, is that they 
predict statistically homogeneous and isotropic Gaussian 
fluctuations with nearly scale invariant power spectrum 

i- 

So far, tests of the scale invariance and/or the Gaus- 
sianity of primordial fluctuations have been done inten- 
sively. In those tests, the statistical homogeneity and 
isotropy are often implicitly assumed. This assump- 
tion, however, should be verifled with observational test. 
For example, the statistical homogeneity of the large- 
scale structure was tested by comparing the means of 
the galaxy distributions in different directions 0]. For 
the CMB temperature held, the anisotropy of the power 
spectrum has rigorously been tested, particularly in the 
context of the hemispherical asymmetry Similar 
tests in real space of the mean, variance, skewness and 
kurtosis were performed in [l^ . 

In this paper we illustrate a method to test the statis- 
tical isotropy of the mean of the CMB temperature fluc- 
tuations in the spherical harmonic space, and we apply it 
to the Wilkinson Microwave Anisotropic Probe (WMAP) 
seven-year observation data. This method can be poten- 
tially useful for the forthcoming PLANCK and future 
CMB surveys, for which contaminations of the CMB 
maps by instrumental noises will take place for higher 
multipole range. 

The CMB temperature field can be expressed as a 



background value Tq with associated fluctuations AT, 

T{h)^To{h) + AT{h) , (1) 

where n is a unit direction vector. Under a null hypoth- 
esis that the background Tq is isotropic, To{n) = Tq and 
the fluctuations AT{n) defined over the full sky can be 
expanded in terms of spherical harmonics as 



with 
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where Ygmin) are the spherical harmonic functions eval- 
uated in the direction n. If the primordial fluctuations 
are statistically homogeneous and isotropic in the mean 
and if they obey the Gaussian distribution, 2^-1-1 spheri- 
cal harmonic coefficients a£mS for each £ are independent 
Gaussian variables with the following properties 



and 
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where (• • •) denotes the ensemble average, Cg is the en- 
semble average power spectrum and 5 is the Kroncckcr 
symbol. 

This null hypothesis that the means of agmS are zero 
[Eq. (j4|)] should be verified with observation. Our aim 
of this paper is to test this null hypothesis against the 
alternative hypothesis that the means of a£„iS arc non- 
zero. One of the difficulties in this test is to decorrelate 
the neighboring modes in harmonic space whose corre- 
lations are caused by a cut sky mask. In order to over- 
come this problem, we use the band-power decorrelation 



2 



analysis following the method in The method is 

to decompose the correlated spectrum to mutually inde- 
pendent band-power spectrum by diagonalizing the co- 
variance matrix. We run Monte Carlo simulations to 
calculate the covariance matrix of the mean distribution, 
including anisotropic instrumental noise and the effect of 
the mask. 

Recently, Armcndariz-Picon tested this zero-mean hy- 
pothesis in a smart analytic way with a symmetric cut sky 
mask to obtain uncorrelated and independent variables, 
with the noise contributions neglected [l^l- The author 
found significant evidence for nonzero means of aimS in a 
particular multipole range. He concluded, however, that 
this evidence is statistically insignificant because the sig- 
nal is high only at one bin among eight different multipole 
bins. The main difference between our analysis and his is 
that we make extensive use of Monte Carlo simulations 
to keep the cosmological information as much as possible 
while taking into account the noise contributions and the 
effect of the WMAP mask. 

This paper is organized as follows. In Sec. II we briefly 
review the main properties of the CMB temperature field, 
its spherical harmonic analysis and the cut sky mask. We 
describe the method of our analysis based on Monte Carlo 
simulations with the WMAP seven-year data in Sec. III. 
Our results are presented in Sec. IV and we provide a 
discussion in Sec. V. We conclude our findings in Sec. 
VI. 



II. TEMPERATURE FIELD AND MASK 

The temperature field ATobs observed with instru- 
ments involves a convolution with the detector beam win- 
dow B and the detector noise ATnoise- In addition, since 
we have to divide the all sky into finite pixels, we intro- 
duce a pixel smoothing kernel K, and write the temper- 
ature field as 



ATobs ^K*[B* AT, 



■cal 



AT, 



noise J 5 



(6) 



where the star denotes the convolution and ATrcai 
is a superposition of CMB temperature field Acmb and 
foreground emission F; 



Ar,eai(ri) = ArcMB(n) + F{n). 



(7) 



The foregrounds consist of, for example, the dust emis- 
sion and synchrotron radiation from our own Galaxy and 
emissions from extragalactic objects. The WMAP team 
provides the sky map data with the foregrounds reduced 
by using appropriate foreground templates [l^. How- 
ever, the cleaning procedure does not completely remove 
the foreground contamination along the galactic disk and 
that from extragalactic objects. To avoid the effect of the 
residual foreground, the contaminated regions have to be 
masked out. The mask is defined by a position-dependent 
weight function M{n) as 



M (n) = 



for contamination region 

1 otherwise. 



(8) 



By definition, M{n)[K*B* F] (n) = 0. Then, by con- 
struction, we get the masked sky ATm which does not 
include foregrounds, 

ATm [n) = M [fi) [K*B* ATcmb + K * ATnoise] (n) ■ 

(9) 

Equations (|H) and ([5]) state that the 2^-1-1 variables a^m 
form a set of normally distributed independent variables 
with a variance Ci. However, in reality, we cannot obtain 
the true aim but only the pseudo aim i<ie^^) from the 
masked sky map. Rewriting Eq. ^ in harmonic space, 
a a^^^ from the masked sky map is given by 



mask 



E 



all sky 



(10) 



where o,^^^^^ is from the all sky CMB map which we can 
never obtain and Mi^rn,t'm' is the convolution matrix of 
the mask. The detector pixel noise is well described by 
a Gaussian distribution with zero mean as shown in the 
WMAP papers [l3,[l^ . Therefore if ensemble averages of 



true aemS are zero, ( a 



all sky 



0, the ensemble averages 



of pseudo aemS are zero, {a'^^^) = 0. The details of 
the convolution matrix are not essential in the present 
analysis, but the matter of importance is that a^mS have 
correlations with neighboring modes. 



III. DATA AND ANALYSIS 
A. Data 

For the main analysis in this paper, we use the WMAP 
seven-year foreground-reduced sky maps [l6j from differ- 
ential assemblies (DAs) Ql, Q2, V1^2, and Wl, W2, 
W3, W4, pixelized in the HEALPix 17| sphere pixeliza- 
tion scheme with the resolution parameter A'side = 512. 
These sky maps are provided at the LAMBDA web site 
[m . We co-added them using inverse noise pixel weight- 
ing and form either individual frequency combined maps 
(Q, V, W) or an overall combined map (Q-t-V+W) to 
increase the signal to noise ratio according to 
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W{i) 



(11) 



where i and j are the pixel and DA indices, respectively. 
The noise weighting function Wj(i) is given by Wj{i) = 
N°^^{i)/aQj, where aoj is the noise root- mean square for 
the jth DA, and N°^^{i) is the number of observations 

of the ith pixel for the jth DA, and W{i) = J^]' 

For the mask, to eliminate foreground contamination, 

we use the extended temperature analysis mask (KQ75y7 

sky mask) [l^ . 

We use a power spectrum Ci file and a beam transfer 

function file to produce realizations of the temperature 

map. The power spectrum Ci we use is the best fitting 
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Ci to the WMAP seven-year data set in the ACDM 
(-|-sz-|-lens) model. The and the beam transfer func- 
tion of each DA are also available at the LAMBDA web 
site. 

Our conclusions are based on the results with the 
overall combined map (Q+V+W). However, we analyze 
also the individual frequency maps (Q, V, W) to exam- 
ine effects of the residual foreground emissions and the 
anisotropic instrumental noise. Because their properties 
are different for frequency bands. 



B. Monte Carlo simulation 

We carry out Monte Carlo simulations to test whether 
the temperature map observed by the WMAP has in- 
deed zero mean in harmonic space. Our analysis consists 
of four steps. First, we prepare 10 000 realizations of 
the sky map from the normally distributed a^mS gener- 
ated from the same underlying power spectrum. Second, 
we apply the cut sky mask to these Gaussian maps and 
the temperature maps obtained by the WMAP, and ex- 
tract a^m"^^ from each map through the harmonic analy- 
sis. Third, for the set of a^^'^^s obtained in the previous 
step, we calculated the means of a^l^^^s for each mul- 
tipole moment £. Finally, we decompose the correlated 
spectrum of the mean to the independent band spectrum 
and calculate p values of the WMAP observed values by 
comparing with the distributions realized by simulations. 

The method of generating Gaussian sky maps consists 
of three steps described in detail below. There are two 
important points in generating Gaussian sky maps. The 
first is that Gaussian maps should be generated using the 
same power spectrum as the WMAP sky map. The sec- 
ond is that we must add random noise taking into account 
the inhomogeneous nature of A'obs to each Gaussian sky 
map. 

Step 1 : Generate the pure Gaussian sky maps 
We prepared 10 000 Gaussian sky maps by HEALPix 
IDL facility isynfast for each DA, which are generated 
from standard normally distributed random aimS from a 
given power spectrum Ct. The Ct inputted is the best- 
fit (ACDM (-hsz-f lens) model ) of the WMAP seven-year 
data set. The resolution of the maps is iVgide = 512 which 
is the same as the resolution of WMAP data we use. We 
use the same random seed array for all DAs. That is, we 
generate random sample of 10 000 different universes. 
Step 2 : Add noise 

According to the WMAP team [l3, [iBl , the instrumen- 
tal noise is expected to exhibit a Gaussian distribution 
and the noise variance in a given pixel is inversely propor- 
tional to the number of observations of that pixel, A'obs- 
The rms noise per observation ctq for each DA is given by 
the WMAP team. Since a Gaussian map generated by 
isynfast contains no noise, the noise should be added on 
it. We add random noise considering the inhomogeneity 
of A'obs to all pure Gaussian maps. 

Step 3 : Combine Gaussian maps 



We combine Gaussian maps whose seed number is the 
same using Eq. (jlip and form either individual frequency 
maps (Q, V, W) or an overall combined map (Q+V+W). 



C. Test Statistic 

Throughout this analysis we deal with real coefficients. 
The real coefiicients aimS and conventional complex co- 
efficients AemS are related to each other by 
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if m > 0. 



We define the mean spectrum as 



Mf 



VWTi 



(12) 



Here Var^ is the average of the variance of aimS calculated 
from the 10 000 realizations. 
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mask,CK 



A/" ^ 2£ + 1 



(13) 



where a'^^^'" represents the multipole coefficient of the 
ath realization and A/" = 10000 is the total number of 
realizations. This mean spectrum is normalized by the 
number of samples in a given £, namely, 2^+1. 

The neighboring i modes of the mean spectra are 
strongly correlated with each other, and the number of 
independent modes is much smaller. To obtain statisti- 
cally independent variables, we define the binned mean 
spectrum Al^ as. 



(14) 



where i is the index of the multipole bin from £ = £^^i to 
£ = £+^i. This binned statistical variable is normalized 

by the sample number in a bin, i.e., X^fcf' (^^ + !)• 
The covariance matrix of A^i is defined by 



K,, 



1 



a— 1 



(15) 



where M.f represents the value of the ith multipole bin 
in the ath realization. K is a real symmetric matrix with 
its dimension equal to the number of multipole bins, N ^ 
and it can be diagonalized by a real unitary matrix [/, 



C/XC/^=diag(Ai, A2,..., AAr) = A, 



(16) 



where A^'s the eigenvalues oi K. A window matrix W is 
defined as till 



(17) 
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where 



(18) 



and 



A-V2^diag(A^^/^A-^/^...,A^^/^). (19) 



We define a new decorrelated statistical variable iSf by 
using the window matrix as 



N 



(20) 



Then the correlation matrix of the variable is diagonal 
and reads 



N 



/2 



where *W denotes the transposed matrix. 



IV. RESULT 



% , (21) 
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FIG. 2. The covariance matrix of the mean spectrum, Mi, 
calculated by Monte Carlo simulations. The central red strip 
is diagonal element array and the surrounding blue strips in- 
dicate negative correlation. Only a part (100 < £ < 200) of 
the covariance matrix whose dimension is 300 is shown here. 



We calculate the mean spectrum, Me, over the multi- 
pole range from i = I to £ = 300 of the WMAP seven- 
year data, which is shown in Fig. [T] Figure [2] shows 
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FIG. 1. The mean spectrum. Me, of the overall combined 
map (Q+V+W), defined in Eq. (fT2)) . 

the covariance matrix of Me calculated for 10 000 real- 
izations made by Monte Carlo simulations. We can see 
that neighboring £ modes are correlated with each other 
up to the range A£ ^ 10. 

Therefore we divide multipole range from £ ~ 1-300 
into the bins whose size is A£ = 6, 10, or 20 to calcu- 
late the binned mean spectrum, Aii. We calculate the 



covariance matrix Kij defined by Eq. psp and the win- 
dow matrix Wij defined by Eq. (|17|) . Then we construct 
decorrelated statistical variables of the WMAP seven- 
year data using Wij . 

The decorrelated band mean spectra of the WMAP 
seven-year data are represented in Fig. [31 In these pan- 
els, the solid line indicates the results from the combined 
map (Q+V-fW), 



N 



(22) 



The box represents the standard deviation for each mul- 
tipole bin, 



N 



/2 



, (23) 



and the bin size are taken as A£ = 6, 10 or 20 in calcu- 
lating the binned mean spectrum, Aii. The color data 
points show the results from the individual frequency 
maps (Q, V, W). 

We can see, in the top panel and the middle panel 
in Fig. [31 that the data points in the multipole range 
£ ^ 210-260 are localized in the negative side as a clus- 
ter despite the fact that these data points are random 
and mutually independent. This negative localization is 
clearly evident in the lower panel in Fig. [31 where the bin 
width is taken as A£ = 20. It is particularly worth noting 



5 



that 12th data point (on the multipole bin i = 221-240) 
indicates the nonzero mean evidence with 99.93% confi- 
dence level [C.L.] {a {p value) = 0.068%) for the com- 
bined map (Q+V+W). Figure [Hand Table |T] show the p 
values of the data points in the lower panel in Fig. [S] 
Clearly, the significance is above at the 99% confidence 
level even for the individual Q, V, W maps, and in fact, 
for the V, W maps the significance becomes even larger 
above the 99.9% C.L. 



V. DISCUSSION 

Since in the WMAP observation data the instrumen- 
tal noise becomes dominant at higher multipoles (see 
Appendix A), We have used the multipole range up to 
£ ~ 300. The method we have described can be useful for 
the forthcoming PLANCKs and future CMB surveys , for 
which contaminations of the CMB maps by the instru- 
mental noises will take place for higher multipole range. 

We can see that the statistics in the above-mentioned 
figures exhibit the same behavior across the three fre- 
quency bands (Q, V, W). Therefore it does not seem 
that these results are due to residual foregrounds and/or 
detector noise characteristics which depends on the fre- 
quency band, but the cosmological CMB signals. 

To eliminate residual foreground emissions more con- 
servatively and confirm these results, we repeat the same 
analysis with a more extensive mask which cuts the re- 
gion in galactic latitude |&| < 30° in addition to KQ75y7 
sky mask [l6j . Again we obtain consistent results, and 
in fact, the significance becomes even larger against the 
zero mean hypothesis as shown in Figs. [S] and [SI The p 
values are summarized in TableU] The p value of the 12th 
multipole bin £ — 221-240 goes down to a = 0.034% from 
0.068%, albeit the standard deviation a becomes larger 
due to the loss of the sky area to analyze. 

Among these analyses of three bin sizes, the results of 
the bin size A£ = 6 keep the original information of the 
CMB signal as much as possible and represent the finest 
structure of the means of aimS. However, we find that 
the most intriguing results are for A£ = 20 because they 
represent the negative localization of the means in the 
multipole range £ = 221-240 most significantly. 

On ground that we have tested n = 15 independent 
bins, we should evaluate significance as a whole to draw 
a conclusion against the null hypothesis. If we would like 
to advocate the anomalousness of the 12th multipole bin 
in the lower panel in Fig. [3] at the significant level atot , 
we should require the a of each individual test to satisfy 



1 - atot = (1 - ay 



(24) 



For 71 = 15 independent bins and atot =5%, this yields 
a = 0.34%, and for atot = 1%, a = 0.067%. Table HI 
represents the values of a of the 12th multipole bin and 
the corresponding values of 1 — atot for the combined map 
(Q+V+W) and the individual Q, V, W maps. Clearly, 
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FIG. 3. The decorrelated band mean spectrum obtained by 
the seven-year WMAP data. The boxes and solid line are 
the result of the overall combined map (Q+V+W) and color 
data points are the results of the individual frequency maps 
(Q, V, W). The horizontal width of boxes indicates the bin 
size A£ = 6 (top panel), 10 (middle panel) and 20 (lower 
panel) for calculating the binned mean spectrum, A4i. 



this evidence against the zero-mean hypothesis is still 
significant, keeping the 99% C.L. (Q+V+W). 

The results we have obtained in this paper would pose 
a challenge to the statistical isotropy and homogene- 
ity hypothesis in cosmology, and they should be exam- 
ined through other independent tests. It will be inter- 
esting to test in real space for scales corresponding to 
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FIG. 4. The p values of the data points in the lower panel in 
Fig. [3 for the bin size A£ = 20. 
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FIG. 5. The decorrelated band mean spectra for two different 
masks. The black solid line (thick boxes) is the result for the 
conventional mask (KQ75y7) and the red dashed line (thin 
boxes) is for the more expansive mask (KQ75y7 + cut |fe| < 
30°). 



TABLE I. The p values for each CMB map with the 
KQ75y7 mask (Q, V, W and Q+V+W) and the extensive 
mask(KQ75y7 + cut |6| < 30°) (Q+V+W only). The bin col- 
umn indicates the multipole bounds for calculation the binned 
mean spectrum, Mi. The p value is anomalously small in 12th 
multipole bin £ = 221-240. 
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44.7% 


181 


200 


48.4% 


45.0% 


50.3% 
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79.2% 
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0.006% 


0.068% 
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68.8% 
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82.0% 
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TABLE II. The values of a of the 12th bin and the corre- 
sponding values of 1 — Qtot. 
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FIG. 6. The p values of data points in Fig. [5] 



220 ^ ( ^ 240. For example, a regional test of the 
power spectrum was performed using a set of circular 
regions in real space [J|, and similar analyses were done 
using different size and shape regions in noi. Rath et al. 
studied scale-dependent non-Gaussianities in a model- 
independent way and found si gni ficant anomalies in the 
muhipole range £ = 120-300 |20l-[2^. Since we arc in- 
trigued with the scales at 220 ^ £ ^ 240, we should 
perform an analysis using finer segmentations. 

Can infiationary cosmology explain the anomalous fea- 
ture found in this paper? Interesting discussions have 
been done, for example, in [23l - [25j . Although infiation 
tends to create homogeneous and isotropic fluctuations 
in the mean and variance, they are on top of the ini- 
tially inhomogeneous and anisotropic background uni- 
verse. Hence the fluctuations could be statistically in- 
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homogeneous and anisotropic in total. Another reason 
would be that our universe is indeed anisotropic: even 
if CMB fluctuations are statistically homogeneous and 
isotropic, the CMB photons travel through the inhomoge- 
neous universe, and the gravitational lensing effects could 
impose statistically anisotropic correlations [l^. In this 
case, the fluctuations arc statistically homogeneous and 
isotropic in ensemble average, and the anomaly we find 
is a merely statistical fluctuation due to a particular re- 
alization of the universe. It is worth studying how prob- 
able the observed inhomogeneous universe can create an 
apparently anisotropic signal in the mean of the CMB 
temperature fluctuations in harmonic space through the 
gravitational lensing effect. 



VI. CONCLUSION 

We have illustrated the method to test the statisti- 
cal isotropy of the mean of the full sky CMB tempera- 
ture fluctuations, which can be potentially useful for the 
forthcoming PLANCK data and future CMB surveys. 

By applying the method to the WMAP seven-year 
observation data, We have found a significant evidence 
for a non-zero mean of spherical harmonic coefficients 
aemS at the multipole bin i = 221-240, for the combined 
Q+Y+W map at the 99.93% confidence level with the 
KQ75y7 mask. Using the more extensive mask, the sig- 
nificance becomes larger, at the 99.96% C.L. As a whole, 
despite the 15 individual tests, this evidence against the 
zero- mean hypothesis is still significant, keeping above 
the 99% C.L. This implies that we should challenge the 
zero-mean assumption and the Universe could contain 
patterns on the characteristic scales. 

However, since there might be unknown systematics, 
independent observations like PLANCK are awaited for 
further understanding. 
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Appendix A: contribution from the instrumental 
noise 

In order to assess the confidence levels of our findings 
we have performed Monte Carlo simulations in which we 
took into account the anisotropic noise and the mask. 
Based on the results of the simulations, we advocated 
anomalies in the means of multipole coefficients aimS. 
Nevertheless you may think that the anomalies are due 
to the anisotropic noise. We cannot rebut completely this 
suspicion because we will never know the exact informa- 
tion of particular realization of the instrumental noise 
contained in the WMAP data. However, it is possible to 
estimate the contribution from the anisotropic noise in 
the mean statistics by performing the same analysis for 
the CMB-free combination maps (W1-W2, VI- V2, Ql- 
Q2) , which are mainly determined by the instrumental 
noise. 

Figure [7] shows the results, over plotted on the lower 
panel in Fig. [3] In order to see the effect of noise for 
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FIG. 7. The result from the CMB-free combination (dashed 
lines), over plotted on the lower panel in Fig. [S] It seems that 
the contribution from the noise is smaller than that of the sig- 
nal in the multipole range we are discussing and the behavior 
of the noise nt £ — 221-240 does not show any anomaly. 

higher multipole range, we extend the range of multipole 
up to ^ = 500. The values indicated by three dashed lines 
are computed as follows. A CMB-free map is defined as 



Tji-j2{p) = Tji{p) - Tj2{p) 



(Al) 



where p and j are the pixel and DA indices respectively. 
For this map Tji-j2{i), we rewrite Eq. ([T4)) as 



M 



-I ^'t 



jl— j2,mask 



/VVar, 



A/ES_,(2^ + i) 



(A2) 



Here Var^ is calculated from the three frequency band 
maps (Q, V, W) by Eq. in order to compare with 
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contribution from the instrumental noise with the signals. 
The noise rms increases by a factor of \pl by executing 
Eq. (jAip because the noises in different DA maps are 
uncorrelated. Therefore Eq. (|A2p is divid by Using 
and the window matrix calculated from the 
individual Q, V, W signal maps, we construct decorre- 
lated variables S]^~^'^ by using Eq. (|2ip . The values of 
Sf^~-''^ represent noise-levels of a single map per DA. 

In Fig. ([7|), we see that the contribution from the 
noise is much smaller than that of the signal in the mul- 
tipole range we are discussing and also the behavior of 
the noise at the multipole bin £ = 221-240 does not show 



any anomaly. Further, in fact, expected noise-level of 
the combined maps should be lower than that of a single 
map per DA by a factor of 1/^2 (Q, V-band) or 1/2 (W- 
band) because two or four maps are combined. Therefore 
we may conclude that our findings arc not attributed to 
the anisotropic instrumental noise, but physically signif- 
icant. 

However, for multipole range £ > 300, clearly the noise 
becomes dominant and the statistics behave indepen- 
dently across three frequency bands (Q, V, W). There- 
fore the results in this paper ware obtained by using the 
multipole range up to £ ~ 300. 
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